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We study a system of electrons interacting through long-range Coulomb forces on a one- 
dimensional lattice, by means of a variational ansatz which is the strong-coupling counterpart of 
the Gutzwiller wave function. Our aim is to describe the quantum analogue of Hubbard's classical 
"generalized Wigner crystal". We first analyse charge ordering in a system of spinless fermions, 
with particular attention to the effects of lattice commensurability. We argue that for a general 
(rational) number of electrons per site n there are three regimes, depending on the relative strength 

V of the long-range Coulomb interaction (as compared to the hopping amplitude t) . For very large 

V the quantum ground state differs little from Hubbard's classical solution, for intermediate to 
large values of V we recover essentially the Wigner crystal of the continuum model, and for small 

V the charge modulation amounts to a small-amplitude charge-density wave. We then include the 
spin degrees of freedom and show that in the Wigner crystal regimes (i.e. for large V^) they are 
coupled by an antiferromagnetic kinetic exchange J, which turns out to be smaller than the energy 
scale governing the charge degrees of freedom. Our results shed new light on the insulating phases 
of organic quasi-lD compounds where the long-range part of the interaction is unscreened, and 
magnetic and charge orderings coexist at low temperatures. 
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I. INTRODUCTION 

Compounds that are poor electrical conductors or in- 
sulators often exhibit charge patterns which differ from 
the homogeneous charge distribution encountered in or- 
dinary metals. The charge ordering often occurs at a rel- 
atively high temperature and is accompanied by a mag- 
netic ordering at a lower temperature. These phenom- 
ena are quite generic, as they are observed in a variety of 
systems. Important examples are the three-dimensional 
nickelates, the layered cuprates at specific commensurate 
fillings, layered molecular crystals such as BEDT-TTF 
radical salts and quasi-one-dimensional systems such as 
the organic TMTTF and DCNQI compounds. From a 
theoretical point of view, the full long-ranged Coulomb 
potential should in principle be included when dealing 
with any interacting insulating phase, due to the inef- 
fectiveness of screening. Accordingly, it is not surprising 
to find wide classes of apparently different compounds 
showing similar behaviours in their insulating phases. In 
this respect, the possibility that the two phenomena of 
charge and spin ordering share a common origin is very 
appealing. Yet, it is not obvious that the two separate 
energy scales can be described within the same physical 
mechanism. 

This raises the question of determining a set of un- 
ambiguous signatures of the long-range Coulomb inter- 
actions, to be compared with the experiments. One such 
signature is of course charge ordering itself, as was first 



recognised by Wigner in the framework of the three- 
dimensional electron gasi In the low-density limit, the 
kinetic energy of the electrons becomes negligible and the 
charges arrange themselves in the form of a Wigner lat- 
tice which minimizes the electrostatic energy. However, it 
is sometimes difficult to determine whether the observed 
charge modulation is induced by long-range forces or by 
an instability of the Fermi surface. Moreover, even if the 
main driving mechanism is given by the Coulomb forces, 
the presence of other effects (elastic, magnetic) in real 
compounds can strongly modify the predicted ordering 
pattern. 

When applying the idea of Wigner crystallization to 
electrons in narrow-band solids, we are faced with the 
additional constraint that the charges must sit on the 
sites of a discrete lattice (the host lattice of atoms, to be 
distinguished from the Wigner crystal of electrons). This 
has been pointed out by Hubbard;^ who considered the 
extreme limit of zero bandwidth, a sensible starting point 
if the dominant energy scale is set by the interactions. 

In the present work, we propose a variational treat- 
ment which allows to go beyond the classical limit by 
including the effects of quantum fluctuations. Our re- 
sults indicate that for small but finite hopping param- 
eters, the charge and spin degrees of freedom are ener- 
getically decoupled: charge crystallization, governed by 
the (large) Coulomb energy, naturally gives rise to anti- 
ferromagnetic spin correlations, characterised by a much 
smaller energy scale. The coexistence of charge and spin 
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ordering at low energies can thus be considered as distinc- 
tive of dominant long-range Coulomb interactions, and 
should therefore be ubiquitous in the insulating phases 
of strongly interacting systems. 

It can be argued that for special fillings, short-range 
models with one or few interaction parameters (e.g. the 
Hubbard model for half-filled bands, or the extended 
Hubbard model for quarter-filled bands) can give a cor- 
rect description of the system. Although they give the 
right answer for the ground state charge configuration 
for special fillings, the strong-coupling analysis presented 
here shows that these models can no longer be applied 
as soon as the filling deviates from such special values, 
since in this case they do not yield the proper classical 
ground state (a system with purely local interactions is 
metallic away from half-filling). Therefore, a complete 
understanding of the behaviour of the insulating phases 
as a function of the density can only be achieved by tak- 
ing into account the long range part of the Coulomb po- 
tential. 

We study a system of fermions in one dimension, inter- 
acting via a local repulsion U and a long-ranged Coulomb 
potential Vm = m being the distance between 

electrons. The corresponding tight -binding Hamiltonian 
on an L-site ring is 



H 



'tT + UD + VW 



(1) 



where we have defined the following dimensionless oper- 
ators, describing respectively electron hopping, on-site 
and long-range repulsion. 
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cf^ {ci^a) creates (destroys) an electron of spin a at site 



ni^a — cl^a^a is the occupation number of this state 
and t is the hopping amplitude between nearest-neighbor 
sites. 

The paper is organized as follows. In the next section 
we introduce our variational wave function, starting from 
Hubbard's solution for the classical limit. In Section ITTll 
we study the phenomenon of charge ordering at different 
commensurate fillings by neglecting the magnetic degrees 
of freedom. These are analysed for a quarter-filled band 
in Section IIVI where contact is also made with exper- 
imental results for organic quasi-one-dimensional com- 
pounds. 



II. VARIATIONAL WAVE FUNCTION FOR 
GENERALIZED WIGNER LATTICES 

The limit studied by Hubbard^ is obtained by set- 
ting the bandwidth to zero and J7 — s- oo. In this case 



the problem is classical (but still non-trivial), and corre- 
sponds to minimising the interaction energy by distribut- 
ing the N electrons as homogeneously as possible over 
the L available lattice sites. The difficulty arises from 
the competition between the different interaction param- 
eters Vm- If the potential is convex, as in the Coulomb 
case, a solution can be found for any commensurate fill- 
ing n = N/L = r/s (r and s integers) in the form of the 
generalized Wigner lattice (GWL), a periodic sequence 
of unit cells of size s, each containing r electrons dis- 
tributed according to a minimum spread of interelectron 
distances. The charge density on site I is either 1 or 
according to the following formula:'^ 



m = [Ir/s] - [{I - l)r/s] 



(5) 



where [. . .] gives the integer part of the argument. 

For finite t quantum fiuctuations will tend to delocal- 
ize the electrons away from the classical configuration. 
Nevertheless, the classical solution is still a good approx- 
imation for small t where most of the charge is located 
on the sites predicted by eq. (jSJ. This argument can be 
made more precise for simple fillings n = 1/s. Moving 
a particle away from the classical configuration costs an 
energy 
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1 



(6) 



A density Sn ~ (t/A)^ will thus be allowed to leak to the 
nearest unoccupied sites. The charge gap is A = 0.39T^ at 
quarter filling (n = 1/2) and is rapidly reduced at lower 
fillings, A ~ 2.4:n^V, since the repulsion that prevents 
the occupation of "wrong" sites weakens as the inter- 
electron distance increases. The range of validity of the 
classical approximation can be estimated by requiring 6n 
to be smaller than a threshold value P. This leads to the 
phenomenological criterion that quantum band motion 
has negligible effects as long as V exceeds the value 



Vgwl 



t 



2AVPn^ 



(7) 



Note that Vgwl ^ n-"^, so that it is harder to stabilize 
a classical GWL at low densities, where Vgwl can be 
very large. Taking for instance P = 0.2, eq. lO gives 
VcwL/t « 8 at n = 1/2 and VcwL/t ~ 500 at n = 1/8^ 
In order to study the ground state of the full Hamilto- 
nian we now use the variational ansatz 



(8) 



introduced previously as the strong coupling counterpart 
of the Gutzwiller wave function in the context of the 
Mott-Hubbard transition^. Here is a dimensionless 
variational parameter, T is the kinetic energy operator 
of eq. (j2Jl and l^'oo) is the ground state for t = 0, i.e. 
Hubbard's classical solution. In the same way as the 
Gutzwiller wave function reduces double occupancy and 
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thus states with high potential energy, the factor e~'^'^ 
suppresses configurations with high kinetic energy. 

It can be demonstrated that the wave function is 
charge ordered and insulating for any finite rj — its charge 
stiffness, or equivalently, its Drude weight, vanishes.^- A 
transition to a metallic state occurs for 77 — s- oo, where 
only the configuration with the lowest kinetic energy, i.e. 
the filled Fermi sea, survives. I^'s) is expected to provide 
a fairly accurate description of charge fluctuations in the 
strong coupling limit, but it can no longer be trusted for 
weak couplings, although it tends to the correct limit as 
the interaction vanishes.^^ 

Our wave function is easy to handle for [/ — s- cx) where 
spin fluctuations are suppressed and the problem be- 
comes equivalent to that of spinless fermions. In this 
case, our variational procedure can be used straightfor- 
wardly to study not only simple fillings such as 71 = 1/2 
but also much more complex GWL ground states such 
as n = 11/47, and it can readily be generalized to higher 
dimensions. For finite values of U, where spin fluctua- 
tions have to be taken into account, the problem becomes 
more complicated. However, in the strong coupling limit 
U,V t the energy scales for charge and spin degrees 
of freedom are very different, and the problem of mag- 
netic ordering can be addressed in terms of a low-energy 
spin Hamiltonian, in fact the simple antiferromagnetic 
Heisenberg model for the case n = 1/2. This will be 
discussed in section Hvl 



III. SPINLESS FERMIONS 

We shall first assume that the on-site repulsion is the 
dominant energy scale, and take the limit U — s- 00. This 
eliminates both double occupancy and mixing of different 
spin configurations. We can therefore suppress the spin 
index in the Hamiltonian (^-Q), and ignore the local 
term (|3J. Transforming to Fourier space leads to the 
following simplified model for spinless electrons: 



(9) 



where the band dispersion relation is = —2tcosk, 
Pq = c^^^Ck is the density fluctuation operator, and 
V{q) = — V^log[2(l — cosg)] (the lattice parameter has 
been set equal to 1). The q sums run over the entire Bril- 
louin zone except q = 0, which ensures charge neutrality. 
To avoid multiple counting of the interaction energy, we 
have taken the limit L 00 with the prescription that 
the range of the interaction potential does not exceed 
half the length of the ringi^ The Coulomb potential has 
a logarithmic singularity at long wavelengths, character- 
istic of the one-dimensional case, and an intrinsic cut-off 
at short wavelengths given by the lattice parameter. 

The wave function ((SJ can be worked out exactly for 
any filling n — r/s. We shall present in full detail the 
formalism for n = 1/2. Results will be given for different 
fillings as well. 



A. Filling n = 1/2 

The classical solution in this case corresponds to alter- 
nating occupied and empty sites. This can be written in 
k-space as: 



1*00) = 



n 

keRBZ 
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J|0) 



(10) 



where the product runs over the reduced Brillouin zone 
|fc| < 7r/2, and |0) is the vacuum for electrons. The 

phase correlator e^"^^ is now diagonal and we can write 
straightforwardly the (normalized) variational wave func- 
tion © as 

\'^b{v))= n ]^(e"''''"4+e''^^4+jio) (11) 

keRBZ 

where the normalization factor is given by N'^ = 
2cosh(277efe) (we have set t = 1 so that energies are now 
expressed in units of the hopping parameter). The best 
variational ground state is obtained by minimising the 
energy functional 



(12) 



with respect to the variational parameter. The kinetic 
energy per unit length is readily evaluated 



{T)/L 



RBZ 



dk , /r, N 

— ek tanh(277efc) 
Ztt 



(13) 



The potential energy can be expressed in terms of the 
structure factor S{q) — (pqp^q) / L, 



{VW)/L=1- I dqS{q)V{q). 
The structure factor has a regular part 



(14) 



dk 



RBZ 



1 + sinh(2?7efc) sinh(2?7efc_g) 
cosh(2?7efe) cosh(2?7efc_g 



(15) 

and a Bragg peak at g = tt corresponding to the periodic 
ordering of the charges: 



^(Tr) = L 



dk 



jiBz 27r cosh(2?7efc 



(16) 



Writing the average density distribution as nil) = n + 
77,cos(7rZ), we see that the intensity of the Bragg peaks in 
the structure factor is proportional to the square of the 
order parameter, S{'k)/L = n^. It is finite for any finite 
?7, so that the wave function \'9b) always represents a 
charge-ordered state. 

In the limit r] ^ one recovers the classical GWL. 
The kinetic energy vanishes and the structure factor be- 
comes S{q) = (l/4)5^_7r, so that equation ifT^ gives 
Ed = V^(7r)/8 = -(F/4) log2. The opposite limit 77 ^ 00 
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yields the Hartree-Fock energy of the liquid phase, which 
is equivalent to treating the interaction to lowest order 
in perturbation theory. Equation l(T^ then gives S{q) = 
\q\/2TT and one readily evaluates Eq{V) — —2/tt — cV, 
with c = 7C(3)/87r2 = 0.10657 (C is the Riemann func- 
tion). For finite rj, one can define the condensation energy 
as the energy gained through charge ordering, Econd = 
Eb — Eq. This quantity is shown in figure ^ as a func- 
tion of V. The variational result (full line) is very close 
to the Hartree-Fock result in the broken-symmetry phase 
(dashed line). The latter is known to be a good approxi- 
mation in the strong coupling regime, and is surprisingly 
accurate in the presence of lattice commensurability, 
which strongly reduces charge fluctuations i& Although no 





-0.02 
-0.04 

-0.06 
-0.08 




1 2 3 4 5 6 
V 



0.5 




FIG. 2: Structure factor in the spinless case for different val- 
ues of the interaction strength at n = 1/2: from top to bot- 
tom, V = 1,2,5, lot. The inset shows the intensity of the 
Bragg peak at tt, proportional to the square of the order pa- 
rameter, versus V. The dashed line is the same quantity eval- 
uated in a model with only nearest neighbour interactions. 
The arrow marks the classical limiting value. 

with refsii, where it is shown that the metallic charac- 
ter of the system is enhanced when going from short (i.e. 
nearest-neighbor) to long-range interactions. 



FIG. 1: Total condensation energy (per unit volume) for spin- 
less fermions at n = 1/2 as a function of the interaction 
strength V (energies are in units of t). 

phase transition occurs within our variational approach, 
there is a crossover between a weak-coupling regime, 
where the density modulation is small as compared to 
the average density, and a strong-coupling (or classical) 
regime, where the order parameter approaches the limit- 
ing value of Hubbard's GWL pattern. Correspondingly, 
the condensation energy is exponentially small at low 
V, behaving like E^ond ^ -0.021/ exp(-27r/y), while 
at large V it can be expanded in powers of t/V as 
Econd « 0.64t - 0M7V - 2.59t^/V + . . .. 

A direct indication on how the charge modulation 
builds up is given by the evolution of the structure factor, 
shown in figureEl AtV — 1 (top line) S{q) is very close to 
the structure factor of non-interacting electrons. Upon 
increasing V, the regular part is progressively reduced, 
and some weight is transferred to the Bragg peak (see the 
inset, full line). Despite a marked onset around V ~ t, 
S{t:) is still far from its classical limit (indicated by an 
arrow) at values as high as V — lOt. The variational 
result for a model with only nearest-neighbour interac- 
tions (dashed line) shows that for equal values of V in 
the charge-ordered regime, the long-range nature of the 
Coulomb potential reduces the amplitude of charge order- 
ing, because the energy cost for creating a local charge 
fluctuation is partly compensated by the interaction pa- 
rameters Vrn at distances m > ln^i This is in agreement 



B. General rational fillings 

For generic commensurate fillings of the form n ~ r/s, 
the s harmonics needed to reproduce the GWL pattern 
produce s Bragg peaks in the Brillouin zone at multiples 
of 27r/s (for even s, the points ±7r of the Brillouin zone are 
equivalent), ft is straightforward to extend the analysis 
presented above for r = 1, s = 2 to larger values of s. 
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FIG. 3: Average occupation of site I for different values of 
the interaction strength at n = 1/8: V = 0.5t (triangles), 
V — 20t (circles) , V — 500t (squares) . Lines are guides to the 
eye. 

Figure |31 shows the variational results for the average 
site occupation at a filling n = 1/8, for different val- 
ues of the interaction strength. Unlike the case n = 1/2, 
which exhibits a single crossover from weak to strong cou- 
pling, three different regimes can now be identified. (1) 
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For V/t = 0.5 (triangles) the modulation of the density 
is weak. (II) For larger values of the couphng strength 
(y/t = 20, circles), the charge modulation is strong, and 
the total wave function is essentially made up of localised 
one-electron wave functions with virtually no overlap, 
like in an ordinary Wigner crystal. The GWL (III), where 
only the sites predicted by Hubbard's solution are occu- 
pied, is achieved only at extremely large V [V/t = 2000, 
squares). 




0.1 1 10 100 1000 




0.1 1 10 100 1000 10000 

V 

FIG. 4: Normalised intensities of the Bragg peaks at multiples 
of qwc = 27rn, versus the interaction strength. From top to 
bottom: n = 1/2, 1/4, 1/8. The full, dashed and dotted lines 
correspond respectively to first, second and third harmonics. 
Note the different scale in the three graphs. 



Figure 01 shows the normalised intensities of the peaks 
in the structure factor at multiples of qwc = 2?™, 
namely S(yqwc)/n'^L, as a function of the coupling 
strength, for n = 1/2 (top), 1/4 (center) and 1/8 (bot- 
tom). We see that at fillings n < 1/2, higher harmonics 
progressively appear in the structure factor when increas- 
ing the coupling, in addition to the main peak at qwc- 
The system goes progressively from the weak coupling 
limit (I), where the modulation can be described in terms 
of a single, exponentially small, Bragg peak, to the GWL 
(III), where all the harmonics in the structure factor sat- 
urate to their classical value. Note that the intermediate 
region (II) becomes broader as n decreases. 



C. From the Wigner lattice to the pinned 
charge— density wave 

By construction, our variational wave function exhibits 
long-range order for any finite positive value of rj (and 
thus for any positive value of V in the case of the long- 
range Coulomb interaction) . The evolution of the charge 
modulation proceeds smoothly from a GWL at very large 

V over an ordinary Wigner crystal at large but not too 
large ^ to a small-amplitude charge-density wave from 
intermediate to small values of V . However, it is obvi- 
ously not clear whether a trial wave function which is 
linked to the exact ground state for t/V = can still be 
trusted in a region where V/t < 1, since in this case it 
would be more appropriate to start from the filled Fermi 
sea. 

To examine this point, we consider the soluble model 
of ID spinless fermions with nearest-neighbor interac- 
tion (parameter V). For an average occupation n = 1/2 
this model is equivalent to the XX Z Heisenberg chain 
in zero magnetic field with exchange constants Jx/ Jz = 
Jy/Jz — 'it/V. For V > 2t (uniaxial regime for the 
spin chain) there is true long-range order, but not for 

V < 2t (planar regime) . The absence of long-range or- 
der for small values of V can also be understood within 
the Fermi gas model, where the (crystal) lattice gener- 
ates an Umklapp term; the latter has to exceed a certain 
critical strength in order to sustain long-range charge 
ordes^. Our wave function predicts long-range order for 

V > Vc ~ 1.3t and a weak first-order transition to a 
metallic phase for V — Vc (see the inset of figure [2J). 
Thus the predictions of our variational ansatz agree qual- 
itatively with the exact phase diagram. 

The question of the existence of true long-range or- 
der for ID lattice fermions interacting through the 1/r 
Coulomb potential, i.e., for the Hamiltonian is more 
subtle and still not completely understood. In the ab- 
sence of a lattice potential, SchulaS has shown, using 
bosonization for small V and an analysis of quantum 
fluctuations about the classical Wigner crystal at large 
V, that the long-distance behavior of the density-density 
correlation function is universally given by 

{pix)piO)) - A cos{Qx)e-^^^°Six) ^ (17) 

where Q = 2kp for the spinless model {4:kp for model 
([T|) with U — 00) and the constant c is independent of 
V. Thus Wigner crystallization in the sense of quasi- 
long-range order exists independently of the interaction 
strength in one dimension. This is in sharp contrast to 
three and two space dimensions, where Wigner crystal- 
lization is known to take place at finite values of the 
interaction strength^'^'^ , namely at ~ 100 and ~ 37, 
respectively {rg is defined as the radius of a sphere con- 
taining one electron in average, expressed in units of the 
Bohr radius, and controls the potential/kinetic energy 
ratio). 

One expects that any additional commensurate lattice 
potential will transform the quasi-long-range order into 
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a true long-range order. This has been confirmed on the 
basis of exact diagonahzationi for large V and n = 1/3. 
For small V, estimates of the effects of Umklapp scatter- 
ing (using bosonization^) are also in line with the exis- 
tence of a finite order parameter, although its value and 
the corresponding lattice pinning energy are expected to 
be extremely small. 

In contrast to bosonization, which cannot provide the 
amplitude of the density-density correlation function [the 
strongly interaction-dependent constant A in Eq. H17|l ]. 
our variational approach allows us to estimate the order 
parameter as a function of both the interaction strength 
V and the average site occupation n. In ID the di- 
mensionless parameter of the continuum theory is rg — 
{a/2n)/aB- Within the tight-binding model 0, the dis- 
persion Cfe close to the band edges can be approximated 
by a parabola with an effective mass m* = /{2ta?). 
Using the unscreened Coulomb potential for V and in- 
serting the band mass m* into the definition of the Bohr 
radius, i.e., as — 2ta/V, one obtains 
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FIG. 5: Phase diagram for the spinless model Q. A charge- 
ordered state persists for all commensurate fillings with three 
characteristic regimes, from top to bottom: (I) generalized 
Wigner lattice, (II) Wigner crystal and (III) small-amplitude 
charge-density wave. 



V 



(18) 



Starting from the Wigner crystal phase (II), where the 
electron wave functions are well separated (cf. fig. |3Jl, 
and varying V, we expect to encounter two crossovers, 
the first to a weakly pinned charge-density wave upon 
decreasing V, the second to a GWL upon increasing V. 
The first crossover towards a weakly modulated ground 
state takes place when the spread a of the single-particle 
wavefunctions equals a given fraction S of the average 
inter-particle distance d — 1/n, namely 



5- 



(19) 



In our lattice problem, the spread of the single-particle 
wave functions can be defined as 



d/2 



(20) 



(the origin has been chosen to coincide with an occupied 
site of the GWL). The Lindemann criterion H19|l is ex- 
pected to locate a broad crossover region rather than a 
true phase transition in the present ID case. Therefore 
the choice of S is somehow arbitrary, and we shall take 
S = 0.25. 

We have calculated the crossover VLind predicted by 
the Lindemann criterion (|19|) for various fillings n — 
1/32, 1/16, 1/8, 1/4 and 1/2. Remarkably, the values ob- 
tained at low fillings {n < 1/4) scale linearly with n, 
as expected from continuum theory in view of eq. (|18|l . 
corresponding to a crossover around « 2.5. The sit- 
uation changes at fillings n = 1/4 and n ~ 1/2, where 
the increasing influence of lattice effects favors localiza- 
tion of the particles, thus reducing the value of Vs at the 
crossover (rs « 2 at n = 1/2, see also refi£). As to the 



crossover towards the GWL, we may use the criterion 
(0, i.e. Vgwl ~ 

The resulting phase diagram is shown in figure|5l where 
the two boundaries represent the two crossover regions 
for rational fillings n = 1/s, s integer. Starting from 
Hubbard's classical GWL at F ^ oo (top) , the inclusion 
of quantum fluctuations at a given density n gives rise 
to two successive crossovers. First, around Vgwl the 
particles begin to spill over to the neighboring lattice 
sites, and the ground state becomes similar to the Wigner 
crystal of the continuum model. Reducing further the 
interaction strength progressively suppresses the higher 
harmonics of the density modulation, so that below the 
second crossover around VLind only a small-amplitude 
^%kp" charge-density wave remains. 

These crossover lines have opposite density depen- 
dences. Vgwl ~ n'~^ decreases with increasing density, 
since reducing the inter-particle distance makes it easier 
to confine particles on a single lattice site. Vj^md ~ n- 
increases with density, since the potential/kinetic energy 
ratio scales as V/{nt) [eq. l(TH|l ]. At n = 1/2, where 
lattice effects are most important, Vgwl ~ VLind and 
the two crossovers become indistinguishable: the Wigner 
crystal phase is wiped out, and the system passes directly 
from the classical GWL to the weakly modulated phase. 
Note that in the spinless case presented here, there is full 
particle-hole symmetry relative to n = 1/2, which allows 
to deduce the phase diagram for 1/2 < n < 1. 



D. Fluctuations of polarization 

An insulator, such as our electronic system with a com- 
mensurate charge-ordered ground state, responds to an 
applied electric field E' by a macroscopic polarization 



7 



P = (X/L), where 



X 



(21) 



is the dipole operator and we have set the electron charge 
equal to unity. The dielectric susceptibility is defined as 
X = P/E in the limit E ^ 0. The upper boundii 

2 



(22) 



links X to the mean-square fluctuation of polarization 
(X^) (we have taken {X) = by symmetry) and to the 
charge gap A (defined here as the minimum excitation 
energy for which the dipole matrix element does not van- 
ish). In a Wigner solid the electrons are located close to 
their classical equilibrium positions and each one con- 
tributes an amount [cf. eq. (|20|l ] to the mean-square 
deviation of polarization, giving 



{X^ 



na 



(23) 



According to the discussion of the previous section, the 
properties of the system in region II (the region of the 
conventional Wigner crystal) should agree with the pre- 
dictions of continuum theory. Therefore, when plotted as 
a function of r^, the quantity (cr/d)^ = n{X'^)/L evalu- 
ated at different fillings should fall on a universal curve. 
The results at n = 1/16 (asterisks), n = 1/8 (circles), 
n = 1/4 (squares) and n = 1/2 (triangles) are plotted in 
fig. showing indeed a scaling behavior up to the value 
VGWL{n), indicated by arrows. On the left side of the 
figure, the crossover between the Wigner crystal phase 
(II) and the weakly modulated regime (I) is signaled by 
a marked upturn of the curve around w 2, which is 
roughly filling independent (except for n — 1/2), and 
agrees with the phcnomcnological Lindcmann criterion 

The scaling behavior breaks down when we enter the 
GWL region (III), i.e. at Vgwl- The onset depends on 
the filling, so that the scaling region is progressively re- 
duced as n increases, and it completely disappears at 
n = 1/2, where lattice effects are most important. Note 
that at this particular filling the fluctuations of polariza- 
tion are systematically reduced with respect to the uni- 
versal curve (the reduction is approximately 30% around 
rg « 2), showing that lattice commensurability favors 
charge localization. Above Vgwl (on the right side of 
the figure), where the charges become localized on indi- 
vidual sites, there is a rapid drop in the fluctuations of 
polarization, and it can be shown following the perturba- 
tive arguments of section ^1 that {a/d)^ ~ 1/(46 r^n^). 



IV. MAGNETIC CORRELATIONS 

The assumption of infinite U implies for the present 
one-dimensional case that the different spin configura- 
tions do not mix and that the magnetic exchange con- 
stant is zero. It follows that from a thermodynamic point 
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FIG. 6: The fluctuations of polarization as a function of 
Vs — V/Ant, for different flUings: n = 1/16 (asterisks), 
n = 1/8 (circles), n = 1/4 (squares) and n = 1/2 (trian- 
gles). Scaling breaks down when entering the GWL phase 
(marked by arrows), which is dominated by lattice effects. 



of view the spins behave as free magnetic moments giv- 
ing rise to a Curie-type magnetic susceptibility. In some 
of the organic materials such a behavior is indeed ob- 
served at temperatures where charge ordering sets in, but 
at lower temperatures antiferromagnetic correlations lead 
to a substantial reduction of the magnetic susceptibility 
and finally to a magnetically ordered ground state. In 
this section we show that our variational approach can 
be readily generalized to finite values of U. We limit 
ourselves to large Coulomb interactions, where the mag- 
netic energy scale turns out to be much smaller than the 
energies associated with the charge ordering, and treat 
explicitly the quarter-filled case (n = 1/2), which is rel- 
evant for some quasi-one-dimcnsional compounds such 
as the TMTTF salts^^ ^j. the DCNQI materials^^ (foj. ^ 

collection of recent results see lefM-) . 

We consider the Hamiltonian for large coupling 
constants, U ^ t,V ^ t. Since the variational ansatz 
(|5|) does not penalize configurations with doubly occu- 
pied sites, we use a refined wave functioniS*i£ 



\^BG)=e-^''e-^^\^, 



(24) 



where l^'oo) is now a linear superposition of states that 
can be visualized in real space as 

{ I C72 C73 ••• ajv )} . (25) 

(ct; is the spin of a localized electron and zeroes stand for 
empty sites). Since the classical ground state is highly 
(2^) degenerate, we do not introduce a priori any special 
magnetic order (e.g. antiferromagnetic or ferromagnetic) 
but we allow for all possible spin configurations. In the 
variational wave function (|24(l the term e~^'^ again con- 
trols the delocalization of charge away from the general- 



ized Wigner lattice (|25|l and 



reduces the weight of 



configurations with doubly occupied sites. As there is no 
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prefered spin configuration in the classical GWL, mag- 
netic correlations will be generated exclusively by quan- 
tum fluctuations. 

The energy of the ground state is obtained by mini- 
mizing the expression 



(26) 



(*cx>|e-''*e-2^^e-''*|*oo) 



with respect to the two variational parameters A and rj. 
Unlike what happens in the spinless closed ana- 

lytical expression for i?BG(A, 77) has not been found. In- 
stead, we use the expansion e"**^ ~ 1 — rjT + ^ij^T^ + . . ., 

which is valid for small r/ (i.e. small t)m^ Each factor T 
produces hopping to neighboring sites. Therefore, in or- 
der to allow for spin exchange, we have to expand at 
least to fourth order in rj. The corresponding hopping 
processes are illustrated in figure [7| They lead to an 



X 



X fl X X 



• X • XXXf|»» 

FIG. 7: Intermediate state contributing to the interaction 
energy at fourth order. Dots/crosses stand for filled/empty 
sites respectively. Notice that the first and third processes 
can give rise to spin exchange. 

exchange of neighboring spins as shown explicitly in Ap- 
pendix A, where the expectation values of the hopping 
term T, of the number of doubly occupied sites D and of 
the long-range part of the interaction W are calculated 
to fourth order in 77. Collecting the results of eqs. l|A7p . 
(|A9() and HA14|I . we obtain the total energy 

Ebg (A, 77) = L [ ~ log 2 + 2rit + (2 log 2 - 1)77^^] 
+ JY. ("^oolS^ ■ S,+2 - \n,n,+2\'^^) , (27) 



where the exchange constant J is given by 

J = -AT^h{l + 3e-^) -Srj'^Ue-'^^ (28) 

-Wt]^ [2 log 2 - 1 + (2 log 2 - ^)] . 

8 

We see that the variational parameters 77 (which con- 
trols the delocalization of charge) and A (which reduces 
double occupancy) can be determined independently of 
the spin configurations contributing to the magnetic 



ground state, which depends only on the sign of the 
exchange constant. For antiferromagnetic exchange the 
ground state energy is given hyi^ 



E 



(tools',; • 5*^+2 - |"-i"i+2 1^*00) = 



-§log2. (29) 



If the parameter rj is small enough (which is certainly 
the case in the Wigner crystal phase), it is essentially 
determined by the minimum of the first term in eq. H27|l . 



77 



(21og2- l)y " 



(30) 



On the other hand, A is chosen such as to maximize the 
exchange constant J, namely 



3t 



277[|C/+(21og2-f )F] 
Inserting these two values into eq. H28() we obtain 



J: 



(21og2-l)2l/2[|[/+(2iog2-f )y] 



The resulting ground state energy per site 

Z^^« - 'T^°g'~ (2iogLi)y - i^"^-' 



(31) 



(32) 



(33) 



consists of three terms, the energy of the classical Wigner 
crystal, the energy gain due to charge delocalization and 
the magnetic exchange energy. 

We see here that the finiteness of U leads to an an- 
tiferromagnetic coupling between the spins, which is a 
direct consequence of the quantum fluctuations in the 
charge-ordered state. Recently it has been shown that 
the spin correlations in the Heisenberg chain decay like 
(— l)('^-')(log \i— — j\m^ln the framework of our 

variational approach, the long-range order in the charge 
sector is therefore accompanied by an algebraic magnetic 
order iS 

The variational analysis presented above for n — 1/2 
indicates that in the region where the approach is valid, 
U ^ t,V ^ t, there is a separation of energy scales 
between the charge and the spin degrees of freedom. 
Recently, charge ordering has been observed in the or- 
ganic chain compounds (TMTTF)2 PFg and (TMTTF)2 
AsFg at Too = and lOOK, respectively. No struc- 
tural transition has been observed so far at Tco and the 
magnetic ordering occurs at much lower temperatures 
(around lOK). It is tempting to associate the charge 
ordering, which has been reported on the basis of both 
NMBiS, and dielectric measurements, with the stabi- 
lization of a Wigner lattice, although for a detailed com- 
parison with experiments structural effects may have to 
be taken into account. Previously, experiments on (DI- 
DCNQI)2Ag have already been interpreted in terms of a 
Wigner crystalii^ 
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It has become popular to describe organic conductors 
in terms of the extended Hubbard model, where only 
on-site and nearest-neighbor Coulomb interactions with 
coupling constants U and Vi , respectively, are taken into 
account.^^'^^'^"^ Informations about the parameter values 
can be obtained by exploiting the well-known sum rule 
relating the optical absorption to the total kinetic en- 
ergy. For (TMTTF)2PF6 exact diagonalization studies 
for the kinetic energy reproduce the observed reduction 
of oscillator strength of 0.73^ (as compared to that of 
non-interacting electrons) using the parameters U ~ 8t 
and Vi — 3t or else U = 6.7t and V = 3M (see Fig. 
4a in refi^). Our variational approach applied to the 
same model would give roughly the same values together 
with an exchange constant Ji = 6t-^/{UVi ) « 0.08t. For 
long-range Coulomb interactions, which suppress the ki- 
netic energy much less, the situation is quite different. In 
fact, for TT, = 1/2 and U = 2V, a. reduction factor of 0.73 
would require a Coulomb interaction strength of y = 9t. 
Nevertheless, the exchange constant, given by eq. I|28|) 
would again be of the order of O.OSt. 

In addition to the simple case presented here, the study 
of more complex fillings can be envisaged. For instance, 
the OWL configuration 110110... at filling n ~ 2/3, 
which is relevant to TTF-TCNQ, intrinsically leads to 
a dimerization of the exchange constants, and could re- 
sult in qualitatively new physics. Our approach may also 
shed new light onto the Akp modulations observed in 
TTF-TCNQ25 and in several 1:2 compounds. Most of 
the structural data exhibit a single harmonic modula- 
tion, so that they could be attributed to region (III) of 
the "phase diagram" of fig. El 



V. CONCLUSION 

We have studied a one-dimensional system of electrons 
interacting through the long-range Coulomb interaction 
on a discrete lattice. Starting from the charge-ordered 
classical limit, we have introduced a variational wave 
function which allows to treat the effects of quantum 
fluctuations. This wave function is always insulating and 
charge ordered at any commensurate filling, and is par- 
ticularly well suited in the limit of strong interactions. 

Although no phase transition arises within our treat- 
ment, the phase diagram can be clearly divided into three 
distinct regions separated by crossovers: (I) the clas- 
sical generalized Wigner lattice in the limit where the 
band motion is negligible, (II) a Wigner crystal regime 
at intermediate couplings and (III) a weakly modulated 
regime at low interaction strengths. Lattice effects ap- 
pear to be particularly important at simple commensu- 
rate fillings, where the charge fluctuations are reduced as 
compared to those of a continuous system. This general 
trend agrees with recent results obtained in two space 
dimensions,^ where it was shown that the critical param- 
eter for Wigner crystallization is dramatically reduced 
from the value Vs ~ 35 in continuous models to « 2 in 



a quarter-filled band. 

Concerning the spin degrees of freedom, exchange pro- 
cesses induced by quantum fluctuations were shown to 
give rise to antiferromagnetic correlations. These develop 
out of the charge ordered configuration, with a much 
lower characteristic energy scale than typical Coulomb 
interaction energies. For the quarter-filled band the 
spin sector has been mapped onto the antiferromagnetic 
Heisenberg chain, which exhibits algebraically decaying 
spin correlations. 

The existence of a charge-ordered ground state, as well 
as the aforementioned decoupling of magnetic/charge 
energy scales, are not unique to long-range interac- 
tions. However, besides quantitative differences between 
the genuinely long-ranged and the ad hoc short-ranged 
models, taking into account the unscreened long-range 
Coulomb potential appears to be crucial in the quest 
of a unified view of the insulator-to-metal transition in 
strongly correlated systems, since it is the only means to 
span all the possible band fillings. In this spirit, the Mott 
transition occurring at the "most" commensurate filling 
n — I and the Wigner crystal at n — ^ are intimately 
connected, and can be viewed as two limiting cases of the 
same general phenomenon. 
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APPENDIX A: VARIATIONAL ENERGY AT 
QUARTER-FILLING 

In this appendix we will detail the computation of 
the ground-state energy for the quarter-filled case n = 
N/L = 1/2, using the trial ground state H24|l . We have 
to evaluate the three terms 

EBciv, A) = -t{f) + U{D) + V{W) , (Al) 

where 

(*oo|e-"*e-2^^e-"*|vl/,,) 

Besides the variational parameters rj and A, the wave 
function l^'oo) itself is not known a priori, but has also to 
be determined by minimizing the energy. It will turn out 
that l^'oo) is the ground state of the ID Heisenberg an- 
tiferromagnet, with spins located on even (or odd) sites. 
We use an expansion in powers of rj and therefore have to 
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assume that V is large enough to keep the electrons close 
to their positions in the classical limit. As the analysis is 
straightforward, we only indicate the main stepsi^i 



1. Kinetic energy 

Expanding to fourth order in rj we obtain for the mean 
value of the kinetic energy 



Y(*oo|r'(e-^''-l)T2|*c 
f 4773((^'^|f2|*^))2 . 



(A3) 



To proceed we have to compute correlation functions 
of the type (*oo|T^™|^'oo) (all contributions with odd 
exponents vanish). The operator T'^™ can be decom- 
posed using Wick's theorem. The expectation value of 
the squared kinetic energy operator (m = 1) is easily 
computed and gives (4'oo|i'^|^oo) = L. The expectation 
value of the fourth power of the kinetic energy is more 
involved, but it is easy to see that only terms with two 
or three contractions contribute. These can be expressed 
in terms of density operators and spin ^ operators Si, 
giving 



i eve7i 



1 

3' 

I even 

(A4) 

The term is cancelled by a corresponding term coming 
from the denominator of eq. ljA2p . so that only "linked 
diagrams" contribute to the expectation value of the ki- 
netic energy. In order to calculate the remaining term in 
eq. (|A3|) , we notice that only the first and third processes 
in fig. [71 contribute. We obtain 

(e-^^-l) T'|*oo) = 

i odd^sa 



(A5) 



I even^sa 



Using again the Wick decomposition and the represen- 
tation of the remaining fermionic operators in terms of 
spin and density operators, we find 

(M'«,|f2(e-^^-l)r2|vI/^) = 

-3(e-^-l) («'oo|45, •5,+2-n,n,+2|^'oo)[A6) 



Collecting all the terms in Eq. (|A3|) we arrive at 

(f) = ~27^L+ff{l+?,e-^) Y (*oo|45,-5,+2-n,n,+2|«'« 

i even 

(A7) 



2. On— site interaction energy 

The expansion in 77 of the mean value of the on-site 
interaction energy gives 

{D) = i(vl/^|f 2i5e-2A6f 2|^^)^4 ^ 0(^6) (A8) 

and can be immediately evaluated by comparison with 
eq. l|X6)l . 

(D) = -^V^e-^^ Y (*oo|4S^. • S',+2 - n,n,+2|*oo) • 

i even 

(A9) 



3. Long— range potential energy 



To compute the expectation value (W) of the long- 
range potential, we first recall that l^'oo) is an eigenstate 
of namely (for an overall neutral system) 



VF|*oo) =-|log2 |*c 



(AlO) 



For W' = W + (L/4)log2 all the terms where W acts 
directly on the wave function are thus suppressed and we 
have to evaluate only four terms to order 77^, 



77 



*oo)}+ 0(776) . (AH) 



In this case all the three processes in fig. [7\ contribute. 
An additional complication is the long-range nature of 
the potential which requires the summation of infinite 
series. The first three terms of the r.h.s. of eq. IjAlip can 
be linked to previously calculated expressions in view of 
the relation 

W^'T|*oo) - (21og2-l)f|*oo). (A12) 
The last term in eq. |IA11|) is found to be 

{^oo\f^e-^^W'e-^^f^\-^oo) - (A13) 
4[(21og2 - 1)L2 + (y - 101og2)L] + 



(f -81og2)e 



-2A 



— 7li71i+2|^'c 



i even 
-2 



Terms proportional to L are again cancelled by terms 
from the denominator. Adding the contribution l|A10l) 
we finally get 



L 



+ [l-21og2+ (y-21og2)e 



-2A 



Y (*oo|4S^, • S'j+2 - 71i7i,+2|*oo) ■ (A14) 
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